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ABSTRACT 


The  effect  of  gravity  levels  up  to  400  times  the  local 
gravity  value  (a/g)  was  experimentally  investigated  in  a 
copper,  rotating,  cylindrical  boiler.   Heat  flux  levels  up 
to  80,000  BTU/hr-ft   were  reached. 

Observed  improvement  in  heat  transfer  at  lov/  heat  flux 
values  with  increasing  gravity  levels  confirmed  the  results 
of  previous  investigations.   At  high  heat  flux  values,  a 
consistent  trend  was  not  observed.   VJith  increasing  gravity 
levels,  some  data  indicated  a  decrease  in  heat  transfer 
while  other  data  indicated  an  increase  in  heat  transfer. 
Changes  in  heat  transfer  with  variation  of  the  liquid 
annulus  thickness  was  noted. 
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TABLE  OF  SYMBOLS 
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A      surface  area  of  evaporator,  ft 
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L 

2 
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g      acceleration  of  gravity,  ft/sec 
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h      heat  transfer  coefficient,  BTU/hr-ft  -°F 
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I .   INTRODUCTION 

The  study  of  the  effect  of  increased  gravity  levels  on 
nucleate  pool  boiling  has  been  underway  for  approximately 
eleven  years.   This  study  has  been  pursued  because  boilers 
and  vapor-generation  devices  are  being  extensively  revised 
because  of  new  requirements  for  their  use  in  spaceflight 
applications  and  in  the  widening  use  of  nuclear  reactors 
as  energy  sources.   These  boilers  are  required  to  operate 
at  high  heat  flux  levels  and  under  variable  gravity  levels. 
The  study  of  the  effects  of  high  gravity  levels  on  the 
boiling  phenomenon  may  also  lead  to  a  new  understanding  of 
the  boiling  heat  transfer  mechanism. 

The  objective  of  this  investigation  was  to  obtain  im- 
proved data  on  nucleate  boiling  of  water  at  high  heat  flux 
levels  and  at  gravity  levels  up  to  400  times  the  local 
gravity  value  (a/g)  .   This  was  done  in  order  to  determine 
if  there  was  a  significant  variation  in  heat  transfer  rates 
in  this  region  with  gravity  level  change  and  if  there  was 
any  reversal  in  the  direction  of  the  variation. 

Costello  and  Tuthill  [1]  studied  the  influence  of  in- 
creased gravity  levels  of  1  to  45  a/g  upon  nucleate  boiling 

heat  transfer  at  heat  flux  levels  between  100,000  and 

2 
200,000  BTU/hr  ft  .   In  this  limited  range  of  study,  in- 
creased gravity  levels  required  larger  values  of  superheat 
for  a  constant  heat  flux.   Heat  transfer  was  decreased  for 
increases  in  acceleration  in  this  region. 


Merte  and  Clark  [2]  studied  the  influence  of  increased 
acceleration,  from  1  to  21  a/g ,  upon  nucleate  boiling  heat 
transfer  in  saturated  water.   Under  pool  boiling  conditions, 
they  detected  an  acceleration  effect  which  was  dependant 

upon  the  heat  flux  level.   For  values  of  heat  flux  less 

2 
than  50,000  BTU/hr  ft  ,  increased  acceleration  decreased 

the  superheat  required.   At  heat  flux  levels  above  50,00  0 

2 
BTU/hr  ft',  increased  acceleration  increased  the  superheat 

required.   They  postulated  two  different  modes  of  heat 

transfer  to  explain  this;  nucleate  boiling  and  natural 

convection  heat  transfer.   At  low  levels  of  heat  flux, 

natural  convection  heat  transfer  governed  heat  removal  and 

increased  gravity  levels  enhanced  heat  transfer  rates.   At 

high  heat  flux  levels,  nucleate  boiling  heat  transfer  was 

predominant  and  increased  gravity  levels  decreased  heat 

transfer  rates. 

In  the  natural  convection  region  Ref.  [3]  recommends 

using  formula 

Nu  =  C(GrPr)m   ,    105  <  (GrPr)  <  2xl07 

for  horizontal  flat  plates  where  C  =  . 54 ,  m  =  . 25 ,  Gr  = 
Grashof  number,  Pr  =  Prandtl  number,  and  Nu  =  Nusselt 
number.   Grashof  number  =  -*■=- — *—   and  is  thus  directly  pro- 
portional  to  the  gravity  level.   This  shows  that  as  gravity 
level  increases,  heat  transfer  rates  in  the  natural  convec- 
tion region  also  increase. 


Adelberg  [4]  compared  gravity  forces  with  three  other 
forces  which  also  influence  the  heat  transfer  rate  in 
nucleate  boiling.   These  forces  were  bubble  dynamic  forces, 
surface  tension,  and  forced  convection  induced  drag.   He 
then  reviewed  existing  literature  for  experimental  data 
which  he  considered  relevant  and  compared  his  analysis 
with  this  data. 

Later,  Adelberg  [5]  developed  three  equations  for  heat 
transfer  in  pool  nucleate  boiling.   The  range  of  applica- 
bility of  each  equation  depended  upon  the  gravity  level. 
From  experimental  data  available,  he  concluded  that  his 
equations  for  low  and  moderate  a/g  levels  agree  favorably 
with  the  limited  data  available.   No  data -was  available  to 
verify  his  equation  for  high  a/g  levels. 

Graham  and  Hendricks  [6]  investigated  the  influence  of 
acceleration  upon  the  mechanism  of  heat  transfer  in  satu- 
rated liquid  boiling  from  1  to  10  a/g.   They  demonstrated 
that  increasing  acceleration  tended  to  shift  the  boiling 
characteristic  curve  in  the  direction  of  improved  heat 
transfer  rates  in  proportion  to  the  1/4  power  of  the  dif- 
ference in  wall  temperature  and  fluid  bulk  temperature. 

Gray,  Marto,  and  Joslyn  [7]  investigated  heat  transfer 

behavior  in  a  rotating  boiler  using  acceleration  up  to  200 

2 

a/g  and  heat  fluxes  up  to  505,000  BTU/hr  ft  .   This  study 

confirmed  that  at  low  heat  fluxes,  increased  acceleration 
increased  heat  transfer  rates;  whereas  at  high  heat  fluxes, 
the  opposite  was  true.   This  was  in  agreement  with  previous 


experimenters.   The  region  where  this  crossover  effect  took 
place  occurred  at  progressively  higher  heat  flux  levels  for 
higher  accelerations.   They  concluded  that  unresolved  ef- 
fects such  as  surface  condition  and  aging  affected  the  data 
at  high  heat  fluxes  as  much  as  did  gravity  variations  from 
1  to  200  a/g. 

Marto  and  Gray  [8]  investigated  high  acceleration  and 
heat  flux  effects  on  nucleate  boiling.  They  achieved  ac- 
celerations up  to  400  a/g  and  heat  fluxes  to  818,000  BTU/hr 

2 

ft  .   Their  results  showed  that  for  the  inception  of  nucleate 

boiling  from  natural  convection,  heat  flux  and  required  wall 
superheat  increased  as  acceleration  increased.   They  also 
verified  that  at  low  heat  fluxes,  increased  acceleration  in- 
creased heat  transfer  rates.   But  at  high  heat  fluxes,  heat 
transfer  rates  were  practically  independent  of  acceleration 
variations.   In  this  region  changes  in  heat  transfer  rates 
were  as  small  in  magnitude  as  the  usual  scatter  in  experi- 
mental data. 

Judd  and  Merte  [9]  studied  the  heat  flux  predictions  of 
six  different  models  proposed  to  explain  the  nucleate  boiling 
phenomenon.   They  showed  that  none  of  the  models  adequately 
predict  nucleate  boiling  heat  flux  for  other  than  saturated 
boiling  conditions  at  standard  gravity.   All  the  models  pre- 
dicted decreasing  values  of  heat  flux  as  acceleration  in- 
creased.  These  observations  suggested  the  omission  of  some 
mechanisms  which  are  inoperative  under  standard  conditions 
of  gravity,  but  very  sensitive  to  the  level  of  subcooling 
and  acceleration. 


They  presented  an  empirical  procedure  for  correlating 
the  results  used  in  their  study.  Their  equation  for  heat 
flux  is 

Q/A  =  0.16  V^3^V/3     (Twall-TB)4/3  (a/g)1/3 


+2.39  x  106  Pr"4-1   ^kU~ 


\r 


g   (e£-ev) 


c 


wall   satv 


L 


hfg 


Metric  units  must  be  used  in  this  equation. 

The  first  term  is  the  contribution  of  natural  convection 

and  the  second  term  the  contribution  of  nucleate  boiling. 

They  showed  that  as  a/q  increases,  (T   ., , -T  )  decreases  when 
J  ^  wall   B 

both  (T   ,,-T   ,  )  and  Q/A  are  constant.   At  low  heat  flux 
wall   satv 

levels,  the  first  term  is  important  since  natural  convection 

dominates.   As  a/g  increases,  the  heat  flux  may  either  in- 

1/3 
crease  or  decrease  depending  upon  the  magnitude  of  (a/g) 

4/3 
compared  to  (T   -|-i~TR)    •   At  high  heat  flux  levels,  the 

first  term  becomes  insignificant  and  the  second  term  domi- 

mates.   In  this  region  Q/A  will  decrease  with  increased  a/g 

since  the  (T   , , -T„ )  term  decreases.   This  suggests  that 
wall   B'  ^  -1 

there  may  be  a  crossover  region  of  heat  flux  where  the  ef- 
fect of  increasing  gravity  level  is  reversed. 
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11 •   EQUIPMENT  CONFIGURATION 

To  conduct  this  high  gravity  nucleate  boiling  investi- 
gation, the  boiler  section  of  an  existing  rotating  thermo- 
syphon  apparatus  was  used. 

Figure  1  is  a  schematic  diagram  of  the  thermosyphon 
assembly  and  associated  equipment.   Figure  2  is  a  cross 
sectional  drawing  of  the  thermosyphon.   Figure  3  is  a 
photograph  of  the  equipment  setup  in  the  laboratory. 

The  equipment  was  the  same  as  that  used  by  Woodard  [10] 
with  the  following  modifications.   The  entire  thermosyphon 
was  manufactured  from  pure  copper  stock  with  the  exception 
of  the  two  end  plates  of  the  evaporator.   These  were  manu- 
factured from  304  stainless  steel.   The  grooves  were  designed 
to  reduce  heat  transfer  from  the  barrel  of  the  evaporator 
to  the  end  plates.   The  difference  in  the  thermal  conduc- 
tivity of  the  two  different  materials  served  the  same  pur- 
pose in  this  design.   The  condenser  end  was  left  open  to  the 
atmosphere  by  removing  the  .25  inch  threaded  plug  in  the  end 
section. 

A.   EVAPORATOR  MANUFACTURE 

The  evaporator  section  was  manufactured  from  pure  copper 
stock  and  the  end  sections  from  304  stainless  steel.   It  had 
an  inside  diameter  of  3.125  inches  and  an  outside  diameter 
of  3.635  inches.   The  length  of  the  copper  section  was  3.125 
inches,  and  when  assembled  the  overall  length  was  3.75  inches 
The  actual  heated  length  was  3.5  inches.   The  heating  coil 
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was  11-gauge  Chromel-A  heater  wire  in  a  0.1875  inch  inconel 
sheath  helically  wound  around  the  evaporator. 

The  procedure  for  bonding  the  stainless  steel  end  plates 
to  the  copper  section  and  bonding  the  heater  wire  to  the 
evaporator  was  an  exacting  process.   Its  purpose  was  to  en- 
sure a  uniform  heat  flux  axially  along  the  heater  surface 
and  to  prevent  electrical  shorting  between  the  heater  wire 
and  other  components.   The  procedure  was  as  follows: 

1.  A  thin  sheet  of  gold  copper  foil  was  tack  welded 
to  the  outside  evaporator  surface  and  each  end 
section. 

2.  The  heater  wire  was  wound  over  the  gold-copper 
foil  with  gold-copper  wire  wound  intermittently 
between  the  coils. 

3.  The  stainless  steel  end  plates  were  then  pressed 
in  place. 

4.  The  entire  assembly  was  brazed  in  a  hydrogen 
atmospherically-controlled  oven . 

5.  Sauereisen  cement  was  air  dried  around  the  out- 
side of  the  heater  wire. 

6.  Asbestos  insulation  was  wrapped  around  the  heater 
wire  and  tied  off  with  glass  tape. 

Power  was  supplied  to  the  heater  coils  by  a  130  volt 
D.C.  source  capable  of  delivering  a  maximum  of  100  amperes. 
This  power  was  transmitted  by  a  brush  assembly  consisting 
of  eight  spring-loaded  graphite  brushes  through  two  bronze 
collector  rings  to  the  heater  wire.   The  bronze  collector 
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rings  were  supported  by  twelve  0.25  inch  stud  bolts. 
Phenolic  spacers  we're  used  to  align  the  rings  and  to  in- 
sulate the  bolts  from  the  rings. 

The  outside  face  of  the  evaporator  was  closed  by  two 
pyrex  glass  view  ports.   The  inner  port  rested  on  a 
neoprene  0-ring  and  was  separated  from  the  outer  port  by 
a  compressed  fiber  gasket.   A  glass  tube  was  fused  to  the 
inner  port  and  passed  through  a  hole  in  the  outer  port. 
This  tube  was  used  for  evacuating  and  filling  the  system 
when  used  as  a  thermosyphon . 

In  order  to  me^^sure  the  temperature  gradient  in  the 
evaporator  wall,  four  holes,  0.06  25  inches  in  diameter  and 
2.375  inches  deep,  were  drilled  radially  in  the  evaporator 
wall.   Four  sheathed  thermocouples  were  inserted  into  these 
holes.   The  center  of  one  hole  was  located  0.0  6  25  inches 
from  the  inner  wall  and  another  hole ' s  center  was  located 
the  same  distance  from  the  outer  wall.   These  holes  were 
located  180  degrees  apart.   The  remaining  two  holes  were 
drilled  so  that  their  centers  were  0.125  inches  from  the 
inner  wall.   These  two  holes  were  also  located  180  degrees 
apart. 

Due  to  the  relative  softness  of  copper,  these  holes  were 
not  drilled  true,  and  it  was  thus  necessary  to  measure  the 
misalignment  of  each  hole.   This  was  done  to  determine  the 
exact  position  of  each  thermocouple  junction.   The  following 
procedure  was  used  on  each  hole  to  determine  its  misalignment 
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A  0.0625  inch  diameter  pin  2.5  inches  in  length  was 
hardened  at  one  end  for  approximately  0.5  inches  to  pre- 
clude any  bending.   This  pin  was  inserted  in  the  drilled 
hole  with  the  hardened  end  outward.   Protruding  from  the 
hole  was  a  0.125  inch  section.   A  precision  square  2.375 
inches  high  was  aligned  to  the  edge  of  the  drilled  hole 
using  the  protruding  pin  end  for  alignment.   When  pro- 
perly aligned,  the  precision  square  was  clamped  in  place 
and  the  metallic  alignment  pin  removed.   Then  a  Tungsten 
heliarc  electrode  approximately  six  inches  in  length  and 
0.0625  inches  in  diameter  was  inserted  into  the  hole.   This 
electrode  was  chosen  because  it  is  manufactured  to  be  pre- 
cision straight.   The  misalignment  between  the  precision 
square  and  the  tungsten  electrode  at  the  top  of  the  block 
was  measured  using  precision  feeler  gages.   Since  the 
electrode  was  assumed  to  be  perfectly  straight,  the  mis- 
alignment of  the  electrode  at  the  top  of  the  precision 
square  was  the  same  as  at  the  bottom  of  the  drilled  hole. 

This  procedure  showed  that  the  ends  of  the  drilled 
holes  were  located  as  follows  from  the  centerline  of  the 
evaporator:  1.610,  1.672,  1.683,  and  1.770  inches,  respec- 
tively.  The  hole  nearest  the  inner  wall  was  only  .016 
inches  from  the  wall.   A  slight  protrusion  was  visible. 
Due  to  the  softness  of  the  copper,  the  drilling  may  have 
caused  the  wall  surface  to  bulge  very  slightly. 
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B.   INSTRUMENTATION 

1 .   Temperature  Measurements 

Temperature  measurements  were  obtained  from  four, 
30  gauge,  copper-constantan  thermocouples.  These  were  in- 
serted in  0.0625  inch  diameter  holes  drilled  2.375  inches 
deep  in  the  evaporator  wall.   They  were  sheathed  in  0.0625 
inches  OD  stainless  steel  tubing. 

The  manufacture  of  these  thermocouples  required  a 
careful  procedure.   First,  the  outer  plastic  insulation 
was  removed  from  the  wires  by  using  a  low  flame  to  melt  the 
insulation.   The  residue  was  then  carefully  scraped  from 
the  inner  wires  using  a  surgical  blade.   Extreme  care  was 
used  so  that  the  insulation  of  the  inner  wires  was  not 
ruptured.   This  outer  insulation  was  removed  from  a  one 
foot  section  of  each  thermocouple. 

Stainless  steel  tubing,  0.063  inches  OD,  was  cut 
into  lengths  of  approximately  five  inches.   The  stripped 
thermocouple  wires  were  inserted  into  the  tubing  and  the 
ends  of  the  wires  were  soft  soldered  to  make  the  junction. 
With  this  junction  made,  the  wires  were  pulled  back  until 
only  approximately  1/4  inch  remained  exposed.   A  coating 
of  Conley-Weld  epoxy  adhesive  was  placed  over  the  exposed 
wires  and  junction.   Then  the  wires  were  again  pulled  back 
until  only  the  epoxy  covered  junction  was  exposed.   The 
epoxy  was  used  to  protect  the  junction  and  to  hold  the 
wires  securely  in  the  tube.   After  the  epoxy  had  set,  fine 
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sandpaper  was  used  to  make  a  smooth  finish  and  reduce  the 
size  of  the  tip  to  that  of  the  tubing. 

The  leads  from  the  four  thermocouples  in  the  evapor- 
ator wall  were  brought  down  the  outside  of  the  condenser 
and  through  the  rear  condenser  flange  to  a  junction  board 
at  the  end  of  the  shaft.   From  the  junction  board,  the  leads 
went  inside  the  main  shaft  to  a  slip  ring  assembly  which 
was  attached  to  the  main  shaft  by  a  coupling. 

The  output  from  the  slip  rings  was  fed  to  a  Hewlett- 
Packard  2010C  Data  Acquisition  System.   This  system  provided 
printed  or  visual  output  of  the  thermocouple  readings. 
2 •   Rotational  Speed  Measurement 

The  rotational  speed  of  the  evaporator  was  measured 
using  a  magnetic  pickup  and  a  Systron  Donner  Counter  1013 
Series.   An  electronic  strobe  light  was  used  as  a  check 

3 .  Condenser  Cooling  Water  Temperature  Measurement 
The  inlet  and  outlet  cooling  water  temperatures  to 

the  condenser  section  were  measured  using  a  Hewlett-Packard 
2801A  Quartz  Thermometer.   The  uncertainty  in  these  readings 
was  ±.1°F. 

4 .  Electric  Power  Measurement 

The  power  into  the  evaporator  section  was  measured 
using  DC  voltage  and  amperage  meters. 

5 .  Cooling  Water  Flow  Rate  Measurement 

The  cooling  water  flow  rate  to  the  condenser  was 
measured  using  a  calibrated  flow  rotameter.   The  calibration 
procedure  for  this  rotameter  is  described  in  Ref.  [11].   The 
uncertainty  in  these  readings  was  ±.5%. 
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6 •   Liquid  Annul us  Thickness  Measurement 

The  liquid  annulus  thickness  was  measured  by  a  paper 
rule  scribed  in  1/16  inch  marks  attached  to  the  viewing  glass 
at  the  evaporator  end.   The  uncertainty  in  these  readings 
was  ±.032  inches. 

C.   OPERATING  PROCEDURE 

The  evaporator  was  filled  with  approximately  200  ml  of 
distilled  water  using  a  syringe.   Filling  was  accomplished 
either  from  the  condenser  end  or  the  evaporator  end  thru  the 
glass  tube.   When  using  the  latter  method,  the  glass  tube 
was  sealed  after  filling.   Cooling  water  flow  was  then  es- 
tablished to  a  30%  reading  on  the  rotameter.   The  drive 
motor  was  then  energized  and  the  desired  RPM  established. 
When  a  liquid  annulus  had  formed,  DC  power  was  applied  to 
the  evaporator  heater.   With  the  desired  RPM  and  power  level 
established,  data  was  recorded  when  thermal  equilibrium  was 
established.   Waiting  time  between  data  recording  ranged 
from  5  to  15  minutes  after  a  change  in  variables.   The  fol- 
lowing data  were  recorded:  RPM,  a/g,  power  supplied  to  the 
heater,  four  thermocouple  readings,  coolant  mass  flow  rate, 
temperature  difference  across  the  condenser,  atmospheric 
pressure,  and  liquid  annulus  thickness. 

Securing  the  machine  was  accomplished  by  first  shutting 
off  power  to  the  evaporator  heater  and  then  shutting  off 
the  drive  motor.   Lastly,  the  cooling  water  was  secured. 
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III.   EXPERIMENTAL  RESULTS 

The  original  method  of  heat  flux  calculation  was  to  use 
the  four  thermocouples  installed  in  the  evaporator  wall  to 
determine  the  radial  temperature  distribution  across  the 
evaporator  wall.   With  this  temperature  distribution  deter- 
mined, the  inside  wall  temperature,  T   ,,,  could  be  computed 

wall 

along  with  the  heat  flux,  Q/A,  as  described  in  Appendix  B. 

Power  level  was  limited  to  6  kilowatts.   Above  this 
level,  steaming  occurred  from  the  condenser  end  and  con- 
siderable vapor  was  lost  to  the  atmosphere.   The  initial 
data  run  revealed  that  with  the  highest  power  level  avail- 
able, the  temperature  variation  across  the  evaporator  was 
too  small  to  be  accurately  determined  by  the  thermocouples. 
This  was  due  to  the  high  thermal  conductivity  of  the  copper 
evaporator,  the  thin  wall,  and  the  limit  on  the  power  level. 
The  errors  in  the  location  of  the  thermocouple  wells  de- 
scribed in  section  II  also  contributed  to  the  uncertainty 
in  the  thermocouple  readings.   Therefore  this  method  of 
calculation  was  abandoned. 

The  alternate  method  used  to  calculate  the  evaporator 
wall  temperature  and  the  heat  flux  was  to  make  a  heat 
balance  on  the  condenser  section.   The  heat  supplied  to 
the  evaporator  equals  the  voltage  times  amperage  supplied 
to  the  heater.   The  heat  to  the  condenser  section  equals 
the  heat  supplied  to  evaporator  minus  convection  and  radia- 
tive heat  losses  to  the  atmosphere.   The  heat  removed  from 
the  condenser  section  equals  the  heat  supplied  from  the 
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evaporator  section  plus  the  heat  generated  by  friction  in 
the  seals  and  bearings  plus  the  heat  generated  by  viscous 
dissipation  in  the  coolant.   These  last  two  sources  of  heat 
generation  can  be  combined  into  the  heat  generated  due  to 
rotation.   An  overall  heat  balance  states  that  the  heat 
supplied  to  the  evaporator  equals  the  heat  removed  from 
the  condenser  minus  the  heat  generated  due  to  rotation. 

The  measurement  of  the  heat  due  to  rotation  was  deter- 
mined at  three  different  RPM's.   Cooling  water  flow  was 
established  at  a  constant  rate  and  electric  power  to  the 
evaporator  heater  was  secured.   When  steady  state -. 
conditions  were  established,  the  temperature  difference  of 
the  cooling  water  and  the  mass  flow  rate  were  recorded. 
The  heat  removed  from  the  condenser  was  calculated  using 
the  mass  flow  rate  and  temperature  difference  across  the 
condenser  as  described  in  Appendix  B. 

Evaporator  wall  temperature,  wall  saturation  temperature, 
and  heat  flux  were  then  calculated  by  the  second  method  de- 
scribed in  Appendix  B.   The  specific  data  taken  were  enumer- 
ated under  "Operating  Procedure".   All  reduced  data  appear 
in  Appendix  E. 

Data  runs  1  thru  5  were  taken  at  constant  RPM  with  vary- 
ing values  of  heat  flux.   Data  taking  began  at  25  G's  up  to 
400  G's.   The  liquid  annulus  thickness  decreased  from  .4 
inches  at  run  1  to  .19  inches  at  run  5.   Data  runs  6  thru  11 
were  also  taken  at  constant  RPM  and  varying  heat  flux.   In 
this  series  run  6  started  at  400  G's  and  decreased  to  25  G's 
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during  run  11.   The  liquid  annulus  thickness  decreased  from 
.62  inches  during  run  6  to  .38  inches  during  run  11.   Data 
runs  12  thru  16  were  taken  at  constant  heat  flux  with  vary- 
ing values  of  RPM.   Liquid  annulus  thickness  decreased  from 
.38  inches  to  .25  inches  during  this  series  of  runs.   Data 

run  17  was  taken  at  a  constant  gravity  level,  25  G's,  and  a 

2 
constant  heat  flux  of  80,000  BTU/hr  ft  .   This  run  was  taken 

to  determine  the  effect  of  liquid  annulus  thickness  of  the 

heat  transfer  coefficient.   Liquid  annulus  thickness  started 

at  .5  inches  and  decreased  to  .125  inches.   Data  run  18  was 

taken  at  400  G's  with  varying  values  of  heat  flux.   Liquid 

annulus  thickness  started  at  .63  inches  and  decreased  to 

.56  inches.   This  run  was  taken  to  determine  if  there  was 

any  effect  upon  the  data  due  to  increasing  the  waiting  time 

after  each  variable  change  to  at  least  15  minutes.   This  was 

well  beyond  the  time  that  steady  state  appeared  to  have  been 

established  from  the  indicator  values. 

Visual  observations  of  the  state  of  the  liquid  in  the 

evaporator  were  made  during  each  data  run.   A  stroboscope 

instrument  was  used  to  light  the  evaporator  interior.   At 

low  gravity  levels  agitated  nucleate  boiling  took  place  in 

the  liquid  annulus  at  all  heat  flux  levels.   The  degree  of 

turbulence  increased  as  heat  flux  increased.   At  high  gravity 

levels  the  nucleate  boiling  process  increased  very  slowly  as 

heat  flux  levels  were  increased.   The  degree  of  turbulence 

was  much  less  at  high  gravity  levels.   Any  effect  of  liquid 
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annulus  thickness  variation  upon  nucleate  boiling  was  not 
readily  apparent  from  visual  observations. 
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IX .   DISCUSSION  OF1  RESULTS 

Figure  4  shows  the  results  of  data  runs  2  thru  5.   RPM 
was  held  constant  and  heat  flux  was  varied.   Gravity  levels 
started  at  50  G's  and  progressed  to  400  G's.   Liquid  annulus 
thickness  was  approximately  .2  inches  during  this  series. 
This  figure  shows  that  at  low  heat  fluxes  heat  transfer  is 
improved  as  gravity  levels  increase.   There  is  no  crossover 
region  discernible,  although  the  data  appears  to  merge  to- 
gether at  high  heat  fluxes. 

Figure  5  shows  the  results  of  data  runs  6  thru  9.   RPM 
was  again  held  constant  and  heat  flux  was  varied.   Liquid 
annulus  thickness  was  approximately  .5  inches.   Data  taking 
started  at  400  G's  and  decreased  to  50  G's.   Again  higher 
gravity  levels  show  improved  heat  transfer  rates  at  low 
values  of  heat  flux.   In  this  figure  there  is  a  definite 
crossover  of  the  100,  200,  and  400  G  curves.   The  data  at 
50  G's  appears  to  be  offset  at  superheats  which  are  too 
high. 

Figure  6  shows  the  results  of  data  runs  11  thru  1G. 
Liquid  annulus  thickness  was  approximately  .38  inches  during 
these  runs.   Heat  flux  was  held  constant  and  RPM  was  varied. 
There  is  a  definite  improvement  in  heat  transfer  with  in- 
creasing gravity  level,  but  no  crossover  region  is  evident. 

Figure  7  shows  the  effect  of  liquid  annulus  change  on 
the  heat  transfer  coefficient.   As  the  thickness  of  the 
annulus  decreases,  the  heat  transfer  coefficient  increases. 
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This  agrees  with  the  results  of  the  investigation  of  MacKenzie 
[12]  and  the  results  of  Ref.  13. 

Figure  8  illustrates  the  effect  of  extending  the  waiting 
time  at  each  data  point  from  5-7  minutes  up  to  15  minutes. 
This  figure  is  a  comparison  between  data  runs  6  and  18.   It 
appears  that  the  longer  the  waiting  time  the  more  reliable 
the  data.   Even  though  the  parameter  indicators  show  that 
steady  state  has  been  established  after  approximately  5 
minutes,  waiting  until  at  least  15  minutes  results  in  a 
larger  temperature  difference  and  a  higher  heat  flux.   Thus 
it  takes  much  longer  than  is  readily  apparent  to  reach 
steady  state  conditions. 

Figure  9  is  a  comparison  between  data  collected  at  25  G's 
during  these  tests  and  data  presented  in  Ref.  2  at  20  G's 
and  in  Ref.  8  at  25  G's.   There  is  good  correlation  between 
the  results  of  Ref.  2,  but  the  data  of  Ref.  8  hiis  higher 
values  of  superheat. 

Figure  10  illustrates  the  data  collected  at  50  G's  with 
that  presented  in  Ref.  8  at  the  same  gravity  level.   The 
50  G  data  from  Ref.  8  is  offset  in  the  direction  of  higher 
superheat.   A  possible  explaination  for  this  is  different 
surface  conditions  and  ageing  effects  between  their  test 
surface  and  the  copper  surface  used  here. 

Figure  11  is  a  comparison  between  the  4  00  G  data  during 
this  test  and  the  400  G  data  presented  in  Ref.  8.   There 
is  good  correlation  between  both  investigations. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

1.  During  nucleate  pool  boiling  at  low  heat  flux 
levels,  there  is  a  definite  improvement  in  heat 
transfer  with  increasing  gravity  levels. 

2.  At  higher  heat  flux  levels  some  data  shows  a 
crossover  trend  while  other  data  appear  to  just 
merge.   No  definite  conclusion  can  be  drawn  from 
the  available  information. 

3.  Increasing  the  waiting  time  after  any  variable 
change  beyond  the  apparent  time  of  steady  state 
enhances  the  reliability  of  the  data. 

4.  Variations  in  the  liquid  annulus  thickness 
cause   changes  in  the  heat  transfer  coefficient. 

B.  RECOMMENDATIONS 

1.  Increase  the  waiting  time  after  any  variable 
change  to  at  least  15  minutes  to  increase  the 
reliability  of  the  data. 

2.  Operate  with  the  machine  closed  to  the  atmo- 
sphere so  that  steaming  is  prevented  and  heat 
flux  levels  can  be  raised  above  those  obtained 
in  this  investigation.   This  will  necessitate 

a  device  for  measuring  the  temperature  or  pres- 
sure in  the  evaporator  vapor  space. 
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3.  Provide  fins  on  the  outside  of  the  condenser  to 
increase  the  heat  transfer  rate  to  the  cooling 
water. 

4.  Increase  the  evaporator  wall  thickness  so  that 
larger  temperature  differences  can  be  realized. 

5.  Install  a  permanent  photographic  installation 
so  that  the  boiling  phenomenon  can  be  visually 
recorded . 
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APPENDIX  A 
CALIBRATION  OF  THERMOCOUPLES 

The  copper-constant  thermocouples  were  calibrated  using 
the  Hewlett-Packard  Data  Acquisition  System,  2010C,  for 
millivolt  readout.   The  following  fixed  point  calibrations 
were  made : 

a.  Ice  point  of  water  (32.0  °F) 

b.  Boiling  point  of  water  (212.0  °F  at  14.696  psia) 

c.  Melting  point  of  tin  (449.4  °F) 

d.  Melting  point  of  lead  (621.1  °F) 

Calibration  of  the  ice  point  of  v/ater  was  accomplished 
by  using  a  Thermos  bottle  filled  with  crushed  ice  and  dis- 
tilled water.   Each  thermocouple  was  inserted  in  this  ice 
bath  and  three  separate  sets  of  data  were  recorded.   The 
average  value  of  these  three  data  sets  was  used. 

For  boiling  point  of  water ,  a  large  test  tube  was  filled 
with  distilled  water,  and  a  loose  fitting  cork  with  a  hole 
in  the  center  was  inserted.   Each  thermocouple  was  inserted 
through  the  hole  in  the  cork  to  approximately  1/2  inch 
above  the  water  level  when  violent  boiling  had  occurred. 
Four  separate  sets  of  data  were  taken  and  the  average  value 
used.   An  uncertainty  for  each  thermocouple  was  also  cal- 
culated using  the  highest  and  lowest  reading.   The  largest 
value  was  ±.009  millivolts.   The  standard  temperature  was 
corrected  for  atmospheric  pressure  variation  using  the 
correction  formula  given  in  Ref.  [7]. 
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Calibration  at  the  melting  point  of  tin  was  accomplished 
by  heating  a  pure  tin  sample  in  a  large  test  tube.   A  small 
diameter  glass  tube  with  one  end  closed  was  inserted  into 
the  sample.   Then  the  thermocouple  to  be  calibrated  was 
placed  in  the  glass  tube.   A  high  range  glass  thermometer 
was  also  inserted  in  the  sample.   The  temperature  of  the 
sample  was  then  raised  above  the  melting  point  and  the  heat 
withdrawn.   The  molten  sample  was  allowed  to  cool  until  the 
freezing  temperature  was  reached.   Then  for  a  discernible 
amount  of  time  the  temperature  stabilized.   At  this  point 
data  was  taken.   This  procedure  was  repeated  six  times 
inorder  to  insure  accurate  results.   The  average  of  the 
six  data  points  was  used  and  an  uncertainty  value  was  then 
calculated  for  each  thermocouple.   The  largest  uncertainty 
value  was  ±.004  millivolts. 

The  same  physical  setup  was  used  to  calibrate  the  thermo- 
couples at  the  melting  point  of  lead  as  was  used  for  tin. 
As  the  lead  was  cooled  in  the  molten  state  to  its  melting 
point,  it  did  not  stablize  for  any  appreciable  time.   There- 
fore a  tape  printout  was  obtained  using  the  minimum  time 
interval  possible.   This  data  was  then  plotted  on  graph 
paper  and  the  area  of  minimum  slope  was  used  as  the  calibra- 
tion point.   The  average  reading  in  this  area  of  minimum 
slope  was  determined  and  uncertainty  value  again  calculated. 
The  largest  uncertainty  value  was  ±.1  millivolts. 
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Reference  14  recommends  a  third  order  equation  of  the 
form 

E  =  at  +  bt2  +  ct3 

for  interpolating  thermocouple  readings  between  100  °C  and 
approximately  300  °C.   Using  this  equation  and  the  calibra- 
tion temperatures,  the  constants  were  determined  using  the 
IBM  360  computer  for  each  thermocouple.   Then  using  this 
result,  a  table  of  thermocouple  readings  versus  degrees 
Centigrade  and  Fahrenheit  was  calculated  for  each  thermo- 
couple using  the  computer. 

The  results  of  the  interpolating  equation  were  plotted 
every  50  °C  inorder  to  compare  these  values  with  the  ex- 
perimental data.   At  the  three  experimental  points  the  un- 
certainty values  were  also  plotted.   The  calculated  results 
compared  well  with  the  experimental  data. 
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APPENDIX  B 
DATA  REDUCTION 

A.   CALCULATION  OF  EVAPORATOR  WALL  SATURATION  PRESSURE 

Evaporator 
Centerline 


liquid  annulus 


//////////////////////////     Wall 


r .  = 


r   = 


h 


Therefore 


evaporator  inner  radius 

distance  from  centerline  to  liquid  annulus 

liquid  annulus  thickness 


r   =  r .  -  h 

a     i 

Squaring  both  sides  of  the  equation 

2      2  2 

r    =  r.   -  2r.h  +  h 
a      i       i 

The  assumptions  in  this  calculation  were  that  solid 
body  rotation  takes  place  and  that  the  density  of  the  two 
phase  mixture  could  be  approximated  by  the  density  of  water 
at  212  °F.   This  type  of  rotation  of  a  fluid  moving  as  a 
solid  about  an  axis  is  called  forced-vortex  motion  [15]. 
The  pressure  field  varies  with  radial  and  axial  position. 
Since  the  machine  was  maintained  in  the  horizontal  position 
there  was  no  axial  variation  in  pressure,  and  since  rota- 
tion was  at  high  RPM,  earth's  gravity  was  neglected.   The 
radial  component  of  the  pressure  becomes: 
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3r      gc 

2 
where   p   =  density  of. liquid  at  212  °F,  Lbm/ft 
L 

w   =  angular  velocity  of  evaporator,  rad/sec 

r   =  radial  position,  ft 

2 

gc  =  conversion  factor,  32.17  f t-Lbm/Lbf-sec 

2 
P   =  pressure,  Lbf/ft 

Integrating  equation  (1)  from  r   to  r .  yields: 

a      i 


AP  =  (r.   -  r   ) 

2gc 

Substituting  for  r    yields: 

a   ■* 


(1) 


2 
p  w 

AP  =  — —  h  (2r.  -  h) 

2gc        1 


or 

2 

p  w 

P   .   =  P  x       +   — —  h  (2r.  ~  h) 
satw    atm    _  i 

2 
where   P   .   =  pressure  at  wall,  Lbf/ft 
satw 

2 

P     =  pressure  in  vapor  space,  Lbf/ft 

B.   CALCULATION  OF  EVAPORATOR  WALL  SATURATION  TEMPERATURE 
The  subroutine  for  the  method  of  least  squares  on  the 
IBM  36  0  computer  was  used  to  find  an  appropriate  polynomial 
to  simulate  the  water  saturation  pressure-temperature  curve 
The  limits  of  14.0  to  50.0  psia  were  used  and  a  fourth 
order  polynomial  gave  excellent  correlation  between  actual 
data  and  a  number  of  selected  test  points. 
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The  polynomial  used  was  as  follows: 

T   .   =  139.435  +  6.926-12P   .   -  .168663P   .  2 
satw  satw  satw 

+  .0024795P   .  3  -  .0000148866P   . 
satw  satw 


where 


T   ,   =  saturation  temperature  at  evaporator  wall,  °F 
satw  ^  l 


C.   CALCULATION  OF  TEMPERATURE  PROFILE 

The  heat  flow  through  the  evaporator  wall  during  steady 
state  conditions  can  be  reduced  to  a  one  dimensional  problem 
if  there  is  no  axial  or  circumferential  heat  conduction. 
Under  these  assumptions  the  steady  state  conduction  equation 
in  cylindrical  coordinates  is 

l   d  dT 

i  ~   (k(T)r  ~)  =  0  (1) 

r   dr        dr 

where   r  =  radial  coordinate 

T  =  temperature 

k  (T)  =  thermal  conductivity  of  copper. 
Integrating  equation  (1)  once  gives 

dT 
k(T)r-^—  =  constant  =  B, 
dr  1 

The  thermal  conductivity  is  assumed  to  vary  linearly 
with  temperature.   Therefore, 

k(T)  =  ko  +  m(T-To)  (2) 
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Data  on  the  conductivity  of  copper  was  obtained  from 
Ref.   16  .   In  the  temperature  range  of  interest,  0  °F  to 
260  °F,  the  conductivity  of  copper  is  a  linear  function  of 
temperature.   This  data  was  plotted  and  ko  and  rn,  the 
slope,  were  determined  graphically. 

ko  =  231.77  BTU/hr-ft-°F 
m   =  -.02222  BTU/hr-f t-°F2 

To  =  0  °F 

Substituting  equation  (2)  into  (1)  and  then  integrating 
again  gives: 

[ko  +  m  (T/2-To)]T  =  B   +  B   lnr 

The  subroutine  for  the  method  of  least  squares  on  the 
IBM  360  computer  was  used  to  find  the  constants  B   and  B  . 
Using  these  constants  and  value  of  the  inner  wall  radius, 
T   , ,  was  extrapolated. 

D.   CALCULATION  OF  HEAT  FLUX,  Q/A 
From  Fourier's  Law: 

Q   =    2tfL    [k(t)r   |£] 

or 


Q   =    2ttL   B, 
1 


The    inside    area   of    the    evaporator    is 


A   =    2-rrr.L 
l 
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Then 


9.  Q    -  ^1  BTU 

A  '  2irr  .L  '  r .   '      hr  ft2 
1      1 


E.   CALCULATION  OF  HEAT  FLUX,  Q/A ,   (SECOND  METHOD) 
The  heat  transfer  rate  out  of  the  condenser  is 

BTU 
Q   .  =  m  CT  (T_  -  T..  )    (gi^-) 
out       L   2     1     hr 

where   m  =  coolant  mass  flow  rate  (r ) 

hr 

pain 

CT  =  coolant  specific  heat  (-rr— — ^=r) 
L  ^  Lbm-°F 

T,  =  coolant  inlet  temperature  (°F) 

T„  =  coolant  exit  temperature  (°F) 

Q   ,  computed  was  then  corrected  for  the  heat  loss  due 

out    L 

to  evaporator  rotation. 

Q     corrected 
Then  Q/A  for  the  evaporator 


Area  of  heated  surface 


Q   ,  corrected 
out 


2-nr  .L 

l 


F.   CALCULATION  OF  EVAPORATOR  WALL  TEMPERATURE,  T   , , 

wall 

(SECOND  METHOD) 
For  steady  state  conditions,  the  heat  flow  equation  is 


2iTkL  (T   -  T   1  .  ) 
q  =  £_ wall_  [1?] 

i     c 
In   — 
r . 

l 
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where   T   =  average  of  the  two  center  thermocouple  reading, 


r   =  radial  distance  to  thermocouples/  ft. 
Then  dividing  by  the  inside  area  of  the  evaporator  gives 


k(T   -  T   _  _  ) 
Q/A  =  c     wal1 

r .  In  — 
1     r . 

1 


Solving  for  T   ,,  yields: 
^       wall 

rc 

ri  (ln  rr]  (Q/A) 

rp        =  t      i 

wall     c 


The  value  of  k  was  assumed  constant  at  226.7,  BTU/hr-ft  °F 
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APPENDIX  C 
UNCERTAINTY  ANALYSIS 

In  order  to  determine  the  reliability  of  the  data  col- 
lected and  to  determine  the  most  probable  error,  an  uncer- 
tainty analysis  was  performed  utilizing  the  Kline  and 
McClintock  technique  [19] .   The  analysis  was  performed  on 
the  data  taken  during  run  1  at  1  KW. 

A.   UNCERTAINTY  IN  THE  CALCULATION  OF  HEAT  FLUX 


^  /-n      ii    m  cp  AT 
Q/A  =  q"  =  I 


r7  „      W*  2    W    2    W.m  2    WA  2 

Wg"  _   ,  m       cp       AT.     ,  A  .    . 


m 


W  =  uncertainty  in  subscripted  quantity 

q"  =  10, 045  BTU/hr-ft2 

cp  =  1.005  BTU/Lbm  °F 

A  =  .2386  ft2 

m  =  30% 

AT  =  2.18  °F 

W  =  ±.5% 
m 

"at  ■  ±a 

W    =  ±.001 
cp 

wA  =  +.00001 

Substituting  into  equation  (1)  yields: 


™Sl   =  .0487  =  4.87' 

q 
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B.   UNCERTAINTY  IN  THE  CALCUALTION  OF  EVAPORATOR  WALL 
TEMPERATURE 

r   In  -?.  q» 

T   ,,  =  T   -  _± li . 

wall     c 


W 
Tw 


=  /  (Wm  )2  +  [?-  (In  —  -  i-  )  w  .] 
/   Tc'     Lk       r.    r.     ri 


„..   r.      2     r.  ,   f£   „ 

vk    r    re      v w,  ; 

c  k2       K 


rc 
r.  In  pj-      2 

+  (— W   )  (2) 

k      q 


r.  =  .130  L  ft 
l 


k  -  226.7  BTU/hr-ft-°F 


Wri  =  ±-001 

r   =  .1398  ft 
c 

W    =  ±2 
wTc 

W,  =  ±.1 

k 

q"  =  10,045  BTU/hr-ft2 


Wrc  =  ±-001 
W  „  =  ±.0487 

q 


Substituting  into  equation  (2)  yields 


w 
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C.   UNCERTAINTY  IN  THE  CALCULATION  OF  WALL  SATURATION 
PRESSURE 

2 

p      =  p     +  P^_  h ( 2  r . -h ) 
satw     atm    2g        1 


+  ^~(2rn_h)Wn)2 
psatw  J        patm       2g     1     P 


+  £|l!  (2r.-h)Ww2)2  +  (f4  M2r.-h)Wgc)2 


.2      .   ...2  2 


+  <f-V    <2%  <2rr2h,wh>  <3> 


r.  =  .1302  ft 

l 

W  =  78.5  rad/sec 

Wgc  =  ±.001 

h  =  .0338  BTU/hr-ft2-°F 

Wp  .   =±1.0 

Wr.  =  ±.001 

g   =32.17  ft/sec2 


w  = 

p 

±.05 

wh  = 

±.032 

L   = 

6  0  Lbm/ft3 

W  o  : 

=  ±.1 

Substituting  into  equation  (3)  yields: 
psatw      3.36 


P   .      2158.12 
satw 


=  .0015  =  .15 


-6 
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D.   UNCERTAINTY  IN  THE  CALCULATION  OF  WALL  SATURATION 
TEMPERATURE 


T   .   =  139.434  +  6.92642  P   ^   -  .168663  P 
satw  satw  satw 

+  .0024795  P   ,  3  -  .0000148866  P   .  4 
satw  satw 


9  T   .  2 

Satw 

WT     4.     =    /     7^  WP     4-  (4) 

•J-satw  9    P  ^satw  v    ' 

satw 

Substituting    into    equation    (4)    yields: 


Jsatw  =   3.35 

T   .      212.95 
satw 

The  largest  uncertainty  appe£irs  in  the  calculation  of 
the  heat  flux.   This  is  due  to  the  error  in  the  mass  flow 
rate  measurement  and  to  the  error  in  the  coolant  flow  tem- 
perature difference  measurement. 
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FIGURE    3 
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FIGURE  4 . 

DATA  RUNS  2  THRU  5 

EFFECT  OF  GRAVITY 
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FIGURE  5. 

DATA  RUNS  6  THRU  9 

EFFECT  OF  GRAVITY 
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FIGURE  6. 
DATA  RUNS  11 
THRU  16 
EFFECT  OF 
GRAVITY 


-10 


10         20 

T      -  T   .    (°F) 
wall     satw 


30 


45 


o 

I 


OJ 


+J 

I 

u 
\ 

D 
En 
CQ 


CM 

lo 


4-> 

c 

Q) 
•H 

U 
•H 
<-M 

m 
a; 
o 
u 

n 

CD 

m 

w 

«j 
u 

E-i 
4-> 

rd 

K 


29 


28 


27 


26 


25 


© — -o 


125 


.25 


375 


Annulus  Thickness  (Inches) 


Gravity  Level,  25  G's       2 
Heat  Flux,  80,000  BTU/hr-ft 


FIGURE  7.   EFFECT  OF  ANNULUS  THICKNESS  ON  HEAT 
TRANSFER  COEFFICIENT 
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FIGURE  8. 

EFFECT  OF  WAITING  TIME 
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FIGURE  9. 

COMPARISON  OF-  2  5  G 
DATA  WITH  20  G 
DATA  [2]  AND  25  G 
DATA  [ 8 ] 
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FIGURE  10.   COMPARISION  OF  50  G  DATA  AND  50  G  DATA  [8] 
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